Dynamical properties of liquid lithium above melting. 
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The dynamical properties of liquid lithium at several thermodynamic states near the triple point 
have been studied within the framework of the mode-coupling theory. We present a self-consistent 
scheme which, starting from the knowledge of the static structural properties of the liquid system, 
allows the theoretical calculation of several single particle and collective dynamical properties. The 
study is complemented by performing Molecular Dynamics simulations and the obtained results are 
compared with the available experimental data. 
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I. INTRODUCTION. 

In the last two decades, considerable progress has been 
achieved, both experimental and theoretical, in the un- 
derstanding of the dynamical properties of liquid metals, 
and in particular the alkali metals because of their sta- 
tus as paradigms of the so called "simple liquids" . On 
the experimental side, and restricting to the liquid alkali 
metals, nearly all elements of this group have. been, stuck 
.ied by neutron scattering experiments: LiEfTj, Nan, RbEl 
and Cstl. Moreover, liquid Li has also been investigate 
,by means of inelastic X-ray scattering experiments!!] I 
Closely linked to this progress, has been the development 
of the computer simulation techniques because of its abil- 
ity to go beyond the scattering experiments, as they allow 
to determine certain time correlation functions which are 
not accesible in a real experiment. Several computer sim- 
ulation studies .have already been performed, mainly for 
alkali metalsEIf££l, although other systems have also been 
studied, e.g., the liquid alkaline-earthsEj or liquid leacfcl 

On the theoretical side, this progress can be linked 
with the development of microscopic theories which pro- 
vide a better understanding of the physical mechanisms 
ocurring behind the dynamics of simple liquidsEH. It was 
quite important to notice that the decay of several time- 
dependent properties can be explained by the interplay 
■of two different dynamical processestTtj. The first one, 
which leads to a rapid initial decay, is due to the effects 
of fast, uncorrelated, short range interactions (collisional 
effects) which can be broadly identified with "binary" col- 
lisions. The second process, which usually leads to a long- 
time tail, can be attributed to the non-linear couplings 
of the dynamical property of interest with slowly vary- 
ing collective variables ("modes") as for instance, den- 
sity fluctuations, currents, etc., and it is referred to as a 
mode-coupling process. Although the resulting analytical 
expressions are rather complicated, it is sometimes possi- 
ble to make some simplifying approximations leading to 
more tractable expressions while retaining the essential 
physical features of the process. In this way, when dealing 
with liquids close to the triple point or supercooled liq- 



uids, the relevant coupling modes are the density fluctu- 
ations, and therefore one can ignore the coupling toj 
rents, leading to much simpler expressions. Sjogren^ 
first applied this theory, with several simplifications, to 
calculate the velocity autocorrelation function and its 
memory function, in liquid argon and rubidium at ther- 
modynamic conditions near their triple points, and ob- 
tained results in qualitative agreement with the corre- 
sponding Molecular Dynamics (MD) computer simula- 
tions. Also, it is worth mentioning the theoretical stud- 
ies carried out by Balucani and coworkersOotil for the 
liquid alkali metals close to their triple points, which 
basically focused on some transport properties, such as 
the self-diffusion coefficient and the shear viscosity, as 
well as other single particle dynamical properties (i.e., 
velocity autocorrelation function, its memory function 
and mean square displacement). From their work, which 
included several simplifications into the mode-coupling 
theory, they asserted that the single particle dynamics 
of the alkali metals is, at least qualitatively, well de- 
scribed by this theory. Based on those ideas, we have 
recently developed a theoretical approaches that allows 
a self-consistent calculation of all the above transport 
and single particle dynamical -properties, and which has 
been applied to liquid lithiumEj and the liquid alkaline- 
earthsta near their triple points, leading to theoretical 
results in reasonable agreement with simulations and ex- 
periment. 

At this point, it must be stressed that within the mode- 
coupling theory, both the single particle and the collec- 
tive dynamical magnitudes are closely interwoven and 
therefore, the application of this theory to any liquid 
system should imply the self-consistent solution of the 
analytical expressions appearing in the theory. However, 
none of the above mentioned theoretical calculations have 
been performed in this way; in fact, the usual practice 
is to take the input dynamical magnitudes, needed for 
the evaluation of the mode-coupling expressions, either 
from MD simulations or from some other theoretical ap- 
proximations. As a clear example of this procedure let 
us consider the intermediate scattering function, F(k,t). 
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This is a central magnitude within the mode-coupling 
theory, which appears in the determination of several 
single particle and collective dynamical properties. How- 
ever, we are not aware of any self-consistent theoretical 
calculation for this magnitude, performed either within 
the mode-coupling theory or some simplified version of 
it. Up to our knowledge, F(k,t) Iias_.always been ob- 
tained either from MD simulationsEJo, or within some 
approximation, in particular, tha viscoelastic model or 
some simplification of it£jll3'EfrE3. 

Furthermore, the application of the self-consistent 
scheme of Ref. [2^ to lithium and the alkaline-earths, 
showed some deficiencies in the calculated magnitudes, 
which were attributed to the inaccuracies of the viscoelas- 
tic intermediate scattering functions, and we suggested 
that theoretical efforts should be directed to obtain a 
more accurate description of F(k,i). 

In this work we present a first step in this direc- 
tion, by proposing a self-consistent determination, within 
the mode-coupling theory, of the intermediate scatter- 
ing function. This is performed through a theoretical 
framewoijk-, which generalizes our previous self-consistent 
approached in order to include both collective and sin- 
gle particle dynamical properties. This new procedure is 
then applied to study the dynamical properties of liquid 
lithium at several thermodynamic states above the triple 
point, far -which thera.are recent experimental inelastic 
neutronErEl and X-rayEl scattering data available. 

The study of liquid Li from a fundamental level poses 
a number of noticeable problems. For instance, the ex- 
perimental determination of its static structure shows 
specific problems related to the small mass of the ion, 
and to-the possible importance of electron-delocalization 
effectscJ. Moreover, the determination of an accurate 
effective interionic pair potential is not easy, and sev- 
eral pseudopotentials have been proposcdES. Theoreti- 
cal calculations, using the Variational Modified Hyper- 
netted Chain (YMHNC) theory of liquidscS, as well as 
MD simulationsEj have shown that a good description 
of the equilibrium properties of liquid lithium can be 
obtained by using an interionic pair potential with no 
adjustable parameters-derived from the Neutral Pseudo 
Atom (NPA) methodEl. 

In all theories for the dynamic structure of liquids, the 
equilibrium static properties are assumed to be known a 
priori. By using the NPA to obtain the interionic pair 
potential, and the VMHNC to obtain the liquid static 
structure, we are able to obtain these input data, needed 
for the calculation of the dynamical properties, from the 
only knowledge of the atomic number of the system and 
its thermodynamic state. 

The paper is organized as follows. In section || we 
describe the theory used for the calculation of the dy- 
namical properties of the system and we propose a self- 
consistent scheme for the evaluation of several single par- 



ticle and collective properties. Section [II gives a brief 
description of the computer simulations and in section |^ 
we present the results obtained when this theory is ap- 
plied to liquid lithium at several thermodynamic states 
above the triple point. Finally we sum up and discuss 
our results. 



II. THEORY. 

A. Collective Dynamics. 

The dynamical properties of liquid lithium have 
been measured by means of inelastic neutron scattering 
(INS)ErH at T=470, 526 and 574 K and also by inelas- 
tic X-ray scattering (LXS)Q'H at T=488 K and T=533 K. 
Whereas the IXS experiment provides information on the 
collective dynamics of the system, i.e., those properties 
described by the dynamic structure factor S(k,oj), the 
INS experiment allows also to obtain information about 
the single particle dynamics, which are enclosed in the 
total dynamic structure factor S , to t(fc,w) 



S to t{k,uj) = 



-S(k, 



-S s (k,u), (1) 



where S s (k, to) is the self dynamic structure factor and <7i 
and <7 C denote the incoherent and coherent neutron scat- 
tering cross sections respectively, which in the present 
case take on the following values: cr;( 7 Li) = 0.68 barns 
and <7 c ( 7 Li) = 0.619 barns. The dynamic structure factor 
may be obtained as 



S(k, u>) = — Re F(k, z — ~iuj) . 

TT 



(2) 



where Re stands for the real part and F(k, z) is the 
Laplace transform of the intermediate scattering func- 
tion, F(k,t), i.e., 



F(k,z) 



dt 



(3) 



Resorting to the_memory function formalism, F (k, z) can 
be expressed asEa 



F(k,z) = S(k) 



n 2 {k) 

z + t{k,z) 



(4) 



where z) is the Laplace transform of the second- 
order memory function T(k,t), and 



ft 2 (k) = k 2 /(3mS(k) 



(5) 



where m is the mass of the particles, [3 is the inverse tem- 
perature times the Boltzmann constant and S(k) is the 
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static structure factor of the liquid. We note that as we 
consider spherically symmetric potentials and homoge- 
nous systems, the dynamical magnitudes to be studied 
here will only depend on the modulus k =| k \. 

In simple liquids, the reason for dealing with memory 
functions, is connected with the development of micro- 
scopic theories where these dynamical magnitudes play a 
key rolaiaEl The time decay of these magnitudes is ex- 
plained by the interplay of two different dynamical pro- 
cesses: one, with a rapid initial decay, due to fast, un- 
corrected short range interactions, and a second process 
(known as a mode-coupling process), with a long-time 
tail, due to the non-linear couplings of the specific dy- 
namical magnitude with slowly varying collective vari- 
ables ("modes") as, for instance, density fluctuations, 
currents, etc. Therefore, according-to these ideas, T(k,t) 
is now descomposed as followsEirO 



r(M) = r B (M) + r M c(M), 



(6) 



where the first term, Ts(fc, t), has a fast time dependence, 
and its initial value and curvature can be related to the 
static structure of the system. The second term, the 
mode-coupling contribution TMc{k,t), aims to take into 
account repeated correlated collisions, starts as t 4 , grows 
to a maximum and then decays rather slowly. We now 
briefly describe both terms and for more details we refer 
the reader to Ref. EstM 



1. The binary-collision term. 

At very short times, the memory function is well de- 
scribed by Ts(k,t) only; moreover, both T(k,t) and 
Ts(k,t) have the same initial value and curvature, which 
are determined Jay the first six moments of the dynamic 
structure factorEl In particular, we have 



T(k,0) = T B (k,0) 



(3m 



+ n 2 + ~f l d (k)-n 2 (k), (?) 



where fio is the Einstein frequency and 7^(fc) stands for 
the second moment of the longitudinal current correla- 
tion function; their respective expressions are 



(8) 



The binary term includes all the contributions to 
T(k,t) to order t 2 . Moreover, as the detailed features of 
the "binary" dynamics of systems with continuous inter- 
atomic potentials are rather poorly known, we resort to 
a semi-phenomcnological approximation that reproduces 
the correct short time expansion. Although several func- 
tions could fulfill this requirement, the most widely used 
in the literature are the squared hyperbolic secant and 
the gaussian. Although both start in the same way, their 
tails are different with the gaussian being somewhat nar- 
rower. Here, for reasons to be discussed below, we use 
the gaussian expression, i.e., 



T B (k,t) =T B (k,0) e 



(10) 



where 77 (fc) is a relaxation time, whose form can be de- 
termined from a short time expansion of the formally 
exact expression of the binary term, and the result is re- 
lated to the sixth moment of S(k,uj). In this way, after 
making the superposition approximation ftp, the three- 
particle distribution function, it is obtainecS 

r B (fc,0) 
rf(fc) 



'i 
3fc 2 
2(3m 



(3m 



3tf + 2 7 < (fc) 



dfe 



-ik-r 



— ik-r~\ 



m 

1 r dk' 

2~pJ W) 



fl(r)(A-V)V(r) + 
[fc a V"VXr)] ff (r) [PV^V>(r)] 



IU A ( r -1 -> 

—k^ik') {[S(k') 1] + [5(| k - k' |) - l] } 



*{7f (fcO^fafc-fc'!)}^. 



(11) 



where summation over repeated indices is implied 
(a, /3, 7 = x, y, z). Here 



lf(k) = ™ / dre-^g(r) V a V^(r) , (12) 



and therefore, the relaxation time can be evaluated from 
the knowledge of the interatomic pair potential and its 
derivatives as well as the static structural functions of 
the liquid system. 



l l M =~f dfe-^gir) (k ■ V)V(r) , (9) 

with ip(r) and g(r) denoting respectively the interatomic 
pair potential and the pair distribution function of the 
liquid system with number density p, and k = k/k. 



2. The mode- coupling component. 

The inclusion of a slowly decaying time tail is a ba- 
sic ingredient in order to obtain, at least, a qualita- 
tive description of T(k,t). A rigorous treatment of this 
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term implies the combined use of kinetic and mode- 
coupling theories. In principle, coupling to several modes 
should be considered (density-density coupling, density- 
longitudinal current coupling and density-transversal 
current coupling) but for thermodynamic conditions near 
the triple point the most important contribution arises 
from the density-density coupling. So, restricting the 
mode-coupling component to the density-density i-cpu- 
pling term only, this contribution can be written asE3 



MC 



(*.*) = /" 



dk> 



pm J (2tt) 6 
k ■ k' c(fc') + k ■ (k - k') c(| k-k'\) 

F{\ k-k> \,t)F(k',t)- F B (\ k-k> \,t)F B (k',t) 

(13) 

where F B (k', t) denotes the binary part of the intermedi- 
ate scattering function, F(k',t). We note that the part 
involving the product of the FB(fc,i)'s has the effect of 
making T(k,t) very small at short times; therefore the 
influence of the approximation made for the F B (k,t) is 
negligible beyond a short time interval, because the in- 
termediate and long time features of the r(fc, t) are con- 
trolled by the terra-involving the product of the F(k, i)'s. 
Following Sj6grenE3, we assume that the ratio between 
F(k, t) and its binary part can be approximated by the 
ratio between their corresponding self-parts, i.e., 



F B {k,t) 



F sB {k,t) 
F s (k,t) 



F(k,t). 



(14) 



where F s (k, t) is the self intermediate scattering function 
and F S B{k,t) stands for its binary part, which is now 
approximated by the ideal gas expression 

F sB (k,t) =F (k,t) =cxp[--i-fc 2 t 2 ]. (15) 

zmp 

Now, once F s (k,t) is specified, the self-consistent eval- 
uation of the above formalism will yield to F(k, t), and, 
through them, to all the collective dynamical properties 
discussed so far. An appropiate expression for F s (k,t) is 
provided by the gaussian approximation 

F.(k,t) = exp[-A) f dr(t-r)Z(r)}, (16) 
rap j 

where Z(t) stands for the normalized velocity autocor- 
relation function (VACF), i.e., Z(t) = (v^v^O)) / (vf) , 
of a tagged particle in the fluid. This approximation 
gives correct results in the limits of both small and large 
wavevectors and it also recovers the correct short time 
behaviour. Therefore, the evaluation of F s (k, t) is closely 
connected with the dynamics of a tagged particle in the 



fluid, which will be calculated as described in the next 
section. 

For completeness we also mention here the viscoelastic 
approximation for F(k, t). This model, which amounts to 
approximate r(fc,i) by an exponentially decaying func- 
tion with a single relaxation timeoEil, has been used in 
previous studies of single particle dynamics within the 
mode-coupling formalism as a simple and convenient ex- 
pression for the intermediate scattering function. In this 
work, we will also use the viscoelastic model, with the 
relaxation time proposed by LoveseyEa, in order to com- 
pare its predictions with those derived from the more 
ellaborated self-consistent formalism proposed above. 



B. Single particle dynamics. 

The self-intermediate scattering function probes the 
single particle dynamics over different scales of length, 
ranging from the hydrodynamic limit (fc — > 0) to the free 
particle limit (fc — > oo). Within the gaussian approxima- 
tion (eq. |l^) , its evaluation requires the previous determi- 
nation of the VACF. This is a central magnitude within 
the single particle dynamics, whose associated transport 
coefficient is the self-diffusion coefficient, D 



D 



1 

(3m Jo 



dt Z{t) 



(17) 



In a similar way as we have proceeded for the interme- 
diate scattering function, F(k,t), the dynamical features 
related to the motion of a tagged particle in the fluid are 
studied by resorting to its first-order memory function, 
K(t), defined as 



Z(t) = 



dt' K(t - t') Z(t') . 



(18) 



Now, i^(tl_ua_ split into binary and mode-coupling 
contributionalJ'Ej, i.e., 



K(t) = K B (t)+K MC {t). 



(19) 



where the first contribution, K B (t), describes single mi- 
correlated binary collisions between a tagged particle and 
another one from the surrounding medium, whereas the 
second one, KMc{t), describes the long-time behaviour 
and originates from long-lasting correlation effects among 
the collisions. Below, we briefly describe both contribu- 
tions. 



1. The binary term. 

At very short times K(t) is well described by K B {t) 
only; moreover both have the same initial value (Oq, the 
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Einstein frecuency) and the same initial time decay (td ) , 
which, within the superposition approximation for the 3- 
body distribution function, is given as 



fig 



P 



, - ., , , df[V a W^(r)]g(r)[V a W^(r)} 

T~J3 oTTl 



dk 1 

6pJ (2?if 



■jfik') [SW-lhfik') . (20) 



This quantity can therefore be evaluated in terms of the 
interatomic pair potential and the static structural func- 
tions of the liquid. A reasonable, scmi-phcnomcnological, 
expression for the binary term is provided by the follow- 
ing gaussian ansatz 



K B (t) = n 2 {) e- t2 ' T i 



(21) 



which allows the calculation of K B (t) from the knowledge 
of the static structural functions only. 



2. The mode- coupling term. 

The presence of a slowly varying part in K (t) is essen- 
tial for the correctrdescription of the dynamics of a tagged 
particle in a fluidtiHl. This is taken into account by the 
mode-coupling contribution, which arises from non-linear 
couplings of the velocity of the particle with other modes 
of the fluid. Although coupling tOpSaMetal modes should 
be considered, it has been shownUEJEJ that for ther- 
modynamic conditions near the triple point, the most 
important contribution arises from the density-density 
coupling and therefore we writer™ 25 ! 



K MC (t) 



P 



24/?7T 3 TO 



dk k 2 c\k) [F s {k,t)F(k,t) 

- F sB {k,t)F B {k,t)) . (22) 

As already mentioned, Balucani et a have stud- 
ied some dynamical properties of the liquid alkali met- 
als by using, basically, the same coupling term as given 
in equation ( |22] ) with the F(k,t) taken either from MD 
simulations or evaluated within the viscoelastic model. 
They found that a substantial contribution to the inte- 
gral in equation ( p2f ) comes from the most slowly de- 
caying density fluctuations, namely, those in the region 
around the first maximum, k p , of S(k) and therefore they 
simplified the evaluation of equation (22) by consider- 
ing only a weighted contribution of k 
gral. This approach was also followed in anothe 



k p to the inte- 



study 



of liquid lithium at several thermodynamic statesc3; now 



by studing the integrand in equation (|22j) they con- 
cluded that besides the k = k p contribution, another 
one from k w 0.65k p should also be taken into account. 
However, in a MD study for liquid sodium, Shimojo et 



have suggested that the whole range of wavevectors 
should be evaluated in order to quantitatively account 
for the mode-coupling effects on the memora functions. 
Therefore, in recent studies of liquid lithiumEj and liquid 
alkaline-earth metalsEj close to the triple point, we have 
incorporated the whole integral in equation ( |22| ) leading 
to an improved description of the single particle proper- 
ties. 



C. Self-consistent procedure. 

By approximating the intermediate scattering func- 
tion, Elk,t), by that obtained within the viscoelastic 
modelEj, we developped a self-consistent schemeEj by 
which we studied the single particle dynamics, as rep- 
resented by the VACF, its memory function, the self- 
diffusion coefficient, and the self dynamic structure fac- 
tor. Now, we extend this scheme to include the self- 
consistent calculation of the collective properties, as 
represented by the intermediate scattering function, its 
second-order memory function and the dynamic struc- 
ture factor. 

We start with an estimation for both the mode- 
coupling component of Kit), (e.g., KMc(t) = 0) and 
the intermediate scattering function, F(k,t) (e.g., the 
viscoelastic model). Using the known values of K B (t) 
and equation (|l^), a VACF total memory function is ob- 
tained which, when taken to equation ([i"8|), gives a nor- 
malized VACF. Now, by using the gaussian approxima- 
tion for F s (k, t), as given by equation (|l6|), the evaluation 
of the integral in equation (53) leads to a new estimate 
for the mode-coupling component of the memory func- 
tion. Finally, this "single particle" loop is iterated until 
self-consistency is achieved between the initial and final 
VACF total memory function, K(t). Now, from the pre- 
viously obtained dynamical magnitudes, and using the 
known values of T B (k,t), we evaluate the mode-coupling 
term, TMc(k, t), as given in equation (|l3|), obtaining a to- 
tal second-order memory function. The resulting values 
are taken to equations (g)- (^) leading to a new estimate 
of the intermediate scattering function, F(k,t). 

Now, with this new value for F(k,t), we return back 
to the "single particle" loop, which is once more it- 
erated towards self-consistency; therefrom, a new esti- 
mate for F(k,t) is obtained and again the whole pro- 
cedure is repeated until self-consistency is achieved for 
both the K(t) and the F(k,t). The practical applica- 
tion of the above scheme has shown that reaching self- 
consistency for the "single particle" loop requires around 
five iterations, whereas the whole process of getting self- 
consistency for both the K{t) and the F(k,t) needs be- 
tween six and ten iterations. 

In this paper, we apply this scheme to study liquid 
lithium at several thermodynamic states. The results ob- 



5 



taincd within this scheme will usually be compared with 
those predicted by the viscoelastic approximation, that 
is, those results obtained using the viscoelastic F(k,t) 
and performing just a "single particle" loop. 

Finally, we end up by signalling that the self-diffusion 
coefficient may be computed either via equation ( [l7|) or 
also from the VACF total memory function, as 



D~ 



m(3 



dt K{t) . 



(23) 



D. Shear viscosity. 

Another interesting transport property is the shear vis- 
cosity coefficient, 77, which can be obtained as the time in- 
tegral of the stress autocorrelation function (SACF), rj(t). 
This latter magnitude stands for the time autocorrelation 
function of the non-diagonal elements of the stress tensor, 
see eqn. ( p^ ) below. Moreover, r](t) can be decomposed 
into three contributions, a purely kinetic term, rjkk(t), a 
purely potential term, r] pp (t), and a crossed term, rjk p {t). 
However, for the liquid range close to the triple point, 
the contributions to, ry coming from the first and last 
terms are negligible^, and therefore in the present calcu- 
lations we just consider r] pp (t) = 7](t), the purely poten- 
tial part of the SACF. Now, this function is again split 
into two contributions, namely, a binary and a mode- 
coupling component respectively, rj(t) — r?s(t) + 7/Mc(i)- 
Again, the binary part is described by means of a gaus- 
sian ansatz, i.e., 



7] B {t) = G P 



(24) 



where G p , the rigidity modulus, is the initial value of 
both rj(t) and wgCt), and t v is their initial time decay. As 

can be computed from 



shown in Ref. 



both G p and r„ 



the knowledge of the interatomic potential and static 
structural functions of the system. The superposition ap- 
proximation for the three-particle distribution function is 
also used in the evaluation of r^. 

The mode-coupling component, rjMc(t), takes again 
into account the coupling with the slowly decaying 
modes, and for the liquid range studied in this work is 
basically dominated by the coupling to density fluctua- 
tions. Its expression, within the mode-coupling formal- 
ism is given byr 3 tr 5 H 3S l, 



1 

60/3tt 2 



dk k 4 



S'(k) 



S 2 (k) 



[F 2 (k,t)~F B (k,t)] , (25) 



where Fg(k,t) is given by equation ( fl4| ) and S'(k) is the 
derivative of the static structure factor with respect to 



k. This mode-coupling integral is evaluated by using the 
F(k,t) and Fsik.t) previously obtained within the self- 
consistent scheme. 



III. MOLECULAR DYNAMICS SIMULATIONS. 

We have simulated liquid 7 Li at the three thermody- 
namic states specified at the Table Q; one state close to 
the triple point (T=470 K) and the other two at slightly 
higher temperatures (T=526 and 574 K). The simula- 
tions have been performed by using 668 (T=470 K), 740 
(T=526 K) and 735 (T=574 K) particles enclosed in a 
cubic box with periodic boundary conditions. For the 
integratioxL.of the equations of motion we used Beeman's 
algorithm^ 4 ], with a time step of 3 fs. The properties have 
been calculated from the configurations generated during 
a run of 10 5 equilibrium time steps after an equilibration 
period of 10 4 time steps. Moreover, in the case of T=470 
and 574 K, we have also evaluated several fc-dependent 
properties; for T=470 K we have considered 20 differ- 
ent fc-values in a fc-range between 0.25 A" 1 and 5 A" 1 , 
whereas for T=574 K we considered 17 different fc-values 
in the same fc-range. 

The computation of both the normalized VACF and 
mean square displacement has been performed according 
to the standard procedures^. The VACF total mem- 
ory fucntion is calculated by solving the equation (fill) 
as follows. From the computed values of the VACF, we 
obtain its time derivative, i.e., the left hand side of equa- 
tion (|l8|). Then, we start with an initial estimation for 
K(t) consisting of 50 points uniformly distributed within 
a range — i max (usually, 1 ps < i max < 2 ps), from which 
a K(t) for all times is constructed by cubic splines. By 
evaluating the right hand side of equation (p^), an esti- 
mation is now obtained for Z(t) which is compared with 
its exact value so as to obtain the mean square deviation 
of the initial estimation for Kit). The software package 
MERLINlj is now used to optimize the 50 values of the 
K(t) by minimizing such a mean square deviation. 

The computation of F(k,t) and F s (k,t) follows from 
its definitionEj, whereas the evaluation of the first and 
second-order memory functions of F(k,t) has been per- 
formed by means of the same approach as previously 
mentioned for the VACF total memory function. The 
dynamic structure factors S(k,u) and S a (k,u}) were ob- 
tained by Fourier transformation of the corresponding 
intermediate scattering functions; however, in the case 
of F s (k,t) for small k, it is not possible to use this 
straightforward procedure, since these functions decay 
very slowly, and another procedure, described in Ref. |Tc| , 
was used to obtain S s (k,ui). 

Finally, the shear viscosity coefficient, 77, has been ob- 
tained by the time integral of the SACF, from the follow- 
ing Green-Kubo-like relation 
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ll = W I ^ <^ 7 (*K 7 (0)> dt = I V (t) , (26) 







where the sum is to be made on the circular permutation 
of the indices aj (xy, yz, zx), V is the volume of the 
system and a ai {t) are the time dependent elements of 
the microscopic stress tensor, namely, 



N N 

a^{t) = m£>f (t)«7(t) + £ r ?i( l )I'7>(t) , (27) 



which is computed in terms of the velocities vf (t), rel- 
ative forces F?At), and relative distances r?-(i) between 
the particles. 



IV. RESULTS. 

We have applied the preceeding theoretical formalism 
to study the dynamical properties of liquid 7 Li at three 
thermodynamic states close to the triple point, which cor- 
respond to those of the neutron scattering experiments 
(see Table I). The input data required by the theory are 
both the interatomic pair potential and its derivatives 
as well as the liquid static structural properties, g(r), 
S'(fc), c{k). We have recently proposed a new interatomic 
pair potential for liquid hthiuna, based on the Neutral 
Pseudo Atom (NPA) methodEa, which, when used in 
conjunction with the Variational Modified Hyperne 
Chain approximation (VMHNC) theory of liquidsLZH 
has proved to be highly accurate for the calculation of 
the liquid static structure and the thermodynamic prop-i 
erties. Moreover, MD simulations using this potential!^! 
have shown that the dynamical properties are also well 
reproduced when compared with the available experi- 
mental data. We will therefore use both the NPA inter- 
atomic pair potential and the corresponding liquid static 
structure results, obtained within the VMHNC theory, 
as the input data required for the calculation of the dy- 
namical properties. We want to stress here that the only 
input parameters needed for the calculation of the inter- 
atomic pair potential and the liquid static structure are 
the atomic number of the system and its thermodynamic 
state (number density and temperature). 



A. Intermediate scattering function. 

The intermediate scattering function, F(k,t), embod- 
ies the information concerning the collective dynamics of 
density fluctuations over both the length and time scales. 
Moreover, it is a basic ingredient for the evaluation of 
both collective and single particle properties; therefore 
it is important to achieve an accurate description of this 
magnitude. 



First, in figures ^ we show, for T=470 K and 574 
K at several fc values, the MD results obtained for the 
first and second-order memory functions of the F(k,t). 
The MD results for the second-order memory function, 
Tmd (&,£), show that for k < k p it remains positive 
at all times, whereas for greater k's it can take nega- 
tive values. It displays a rapidly decaying part at small 
times and a long-time tail which becomes longer for the 
smaller fc-values, and it exhibits an oscillatory behaviour 
for k > 1.3 w kp/2. On the other hand, the correspond- 
ing MD results for the first-order memory function, ex- 
hibit a negative minimum at all the fc-values considered, 
although its absolute value decreases for increasing fc's 
and becomes rather small for k fc p ; moreover it ex- 
hibits weak oscillations for the small fc-values. In fact, 
very similar Qualitative features were also obtained by 
Shimojo et a&j in their MD study for liquid Na near the 
triple point. 

Figures also show the theoretical T(k,t) obtai ned 
within the iterative scheme described in section II C . As 



already mentioned, a number of around eight iterations 
were enough to achieve self-consistency for this magni- 
tude. First, we note that the theoretical T(fc,i) repro- 
duce, qualitatively, the corresponding MD results. In 
particular, the short time behaviour, which is dominated 
by the binary component, is very well described. This 
is not unexpected, as the present theoretical formalism 
imposes the exact initial values of both r(fc, t) and its 
second derivative through equations (Q) and ( p"i| ) respec- 
tively. Moreover, it is important to note that the overall 
amplitude of the decaying tail of T(fc,i) is also, in gen- 
eral, well described by the present formalism. The main 
discrepancies with the MD results are the amplitudes of 
the oscillations around the decaying tail, which occur for 
intermediate times where the mode-coupling component, 
TMc(k,t), dominates. We consider that the main rea- 
son relies on the expression taken for the TMc(k,t), see 
equation ([l3]), where we have only included the density- 
density coupling term. Although with this term only, it is 
possible to reproduce the main features of the Tmd^, t), 
further improvements would in principle require the in- 
clusion of other terms, especially those related to the 
density-currents couplings. Furthermore, given the self- 
consistent nature of the present theoretical framework, 
the approximations made in the description of other mag- 
nitudes, such as F s (k,t) and Z(t), will also exert some 
influence and therefore, a more rigorous treatment of 
these magnitudes would improve the theoretical results 
obtained for the r(fc, t). 

On the other hand, we note that the gaussian ansatz 
chosen for the description of the binary part, leads to 
theoretical results in reasonable agreement with the sim- 
ulation results. We have verified that the choice of the 
squared hyperbolic secant in equation ( |l0|) leads to a 
worse description of the short time behaviour of the 
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r(fc, t), with the binary term becoming clearly wider than 
the MD results for the large fc values. 

In figures we have also included., for comparison, 
the results of the viscoelastic modelE] for the second- 
order memory function, T v i sc (k,t). In this model, this 
function is represented by a simple exponential decaying 
function where the only exact value is the initial value; 
this explains the strong differences found between the 
F v is C (k,t) and the corresponding MD results. However, 
it is interesting to observe that for small fc-values the vis- 
coelastic model underestimates the MD data at all times, 
whereas around the first peak and for higher wavevec- 
tors there is a compensation between the short time be- 
haviour, where the viscoelastic model is too low, and the 
intermediate times, where it is too high when compared 
with the MD results. 

Figures §]| show, for T=470 K and 574 K at several k 
values, the corresponding MD results for the intermediate 
scattering functions, FMn{k, i). It is observed that the 
FMT>{k, t) exhibit an oscillatory behaviour which persists 
until k w 2 with the amplitude of the oscillations 

being stronger for the smaller fc-values. These oscillations 
do not take place around zero, but instead the behaviour 
of Fy\x)(k,t) after the first minimum looks like an oscil- 
latory one superimposed on a monotonically decreasing 
tail. 

These characteristics are qualitatively reproduced by 
the theoretical F(k,t) obtained within the present for- 
malism. In particular, the amplitude of the decaying tail 
is rather well reproduced, while the main discrepancies 
are located for the small fc-values where the amplitude 
of the oscillations is underestimated. This is an impor- 
tant improvements-over the results obtained within the 
viscoelastic modeE^, since the oscillations of F v ; sc (fc,i) 
take place basically around zero. This has important 
consequences on the single-particle dynamical properties 
at intermediate times, as will be shown in the next sec- 
tion. The deficiencies of F v i sc (k,t) are, however, not so 
marked as those previously seen for the corresponding 
second-order memory function. This improvement is ex- 
plained because the viscoelastic model incorporates the 
exact initial values of F(k,t), and its second and fourth 
derivatives. Comparison with the MD results shows that 
for the small fc-values, i.e., fc rj 0.25 A -1 , the oscillations 
in the F v i sc (k,t) are strongly damped, for the intermedi- 
ate fc-values, i.e., fc w 1 A -1 and fc rj 1.7 A -1 , the os- 
cillations are stronger and they persist for longer times, 
whereas for those fc-values around fc p the results obtained 
by the viscoelastic model show a good agreement with 
the MD results. This agreement for the regions around 
the main peak can be understood in terms of the com- 
pensation between the short-time and intermediate-time 
deficiencies of the viscoelastic second order memory func- 
tion. _By taking into account the suggestions of Balucani 
et a&3 already mentioned in section II B, this good de- 



scription of F(k,i) provided by the viscoelastic model in 
the region of the main peak, expl aieg why those single 
particle dynamics calculationi3c3~E3 based on the vis- 
coelastic model lead to results in reasonable agreement 
with the corresponding MD and experimental results. 

We stress that, to our knowledge, this is the first the- 
oretical study where the intermediate scattering func- 
tion, F(k, t), has been evaluated directly from its memory 
function/mode-coupling formalism, without resorting to 
any parameters or fitting to any assumed shape. 



B. Self diffusion. 

Another magnitude obtained within the present self- 
consistent scheme is the VACF total memory function, 
and in figure [s| we show, for the three thermodynamic 
states studied in this work, the obtained theoretical total 
memory function, K(t), along with the mode-coupling 
component. Whereas the binary part dominates the be- 
haviour of K {t) for t < 0.06 ps, for longer times the mode- 
coupling part completely determines the shape of K(t). 
For comparison, we have also plotted the total memory 
function obtained from the viscoelastic approximation, 
-R'visc(i)) as well as the MD memory function, KMu(t), 
obtained when the MD normalized VACF is used as in- 



put data in equation (18). First, we note that for T=470 



K the present theoretical K{t) shows a rather good agree- 
ment with the MD results, especially at the longer times 
where the wiggles and the tail of the memory function 
are well reproduced. This represents a significant im- 
provement over previous resultsc3, where the viscoelastic 
approximation was used and, as shown in figure ^, the 
long-lasting tail was not qualitatively reproduced by the 
K V i SC {t). In fact, the improvement of the present theoret- 
ical approach over the viscoelastic approximation can be 
traced back to the different intermediate scattering func- 
tions, F(Jc,t), obtained within both theoretical schemes. 
By analizing which wavevectors are really relevant in the 
integral appearing in equation ( p2| ) it is found that, for 
any time, the dominant contribution comes from two sets 
of fc-values centered around fci rj 1.7 A -1 and k 2 ~ k p , 
where k p (ss 2.5 A -1 ) is the position of the main peak of 
the S(k). By comparing, for these two sets of fc-values, 
the different terms appearing in the integrand with its 
MD counterparts we find that the main discrepancies 
come from the intermediate scattering function, F(k,t) 
at fc « fci . As shown in figure |[ both the present theoret- 
ical and the MD-based F(k S3 k\,t) show an oscillatory 
behaviour for t < 0.20 ps and then decrease towards zero 
taking positive values, whereas F v i sc (k w k\,t) keeps an 
oscillatory behaviour around zero for longer times and 
this is the main reason for the discrepancies among the 
corresponding memory functions observed in figure o for 
t > 0.20 ps. In particular, the deep minimum exhibited 



by the K v i sc (t) at t « 0.2 ps is connected to the fact that 
the F v i sc (k w fci, f) shows at around t w 0.2 ps a negative 
oscillation. 

To examine more closely the adequacy of the mode- 
coupling expression used in equation (p2j), we have also 
evaluated this expression, for T=470 K and 574 K, by 
using the MD results for both F(k,t) and F s (k,t). The 
obtained results, denoted as K^?\t), are shown in the 
insets of figure |5j, where we have also included the cor- 
responding i^MD^)- For t < 0.08 ps, a direct compar- 
ison between them is not possible due to the contribu- 
tion of the binary term to the total memory function. 
For t > 0.08 ps, where the binary part is almost cero, 
we obtain that in the case of T=470 K, the K^\t) 
follows the behaviour of KMv(t) although it is slightly 
higher; this is in qualitative agreement with the sugges- 
tion that, near the triple point, the effect of the other 
coupling terms not included in equation (22), gives rise 
to a small, although negative, contribution to the total 
memory function. On the other hand, the comparison 
of the K^? ] (t) with both K MC (t) and K™g(t) reveals 
that the present theoretical framework leads to an im- 
portant improvement over the results derived from the 
viscoelastic approximation, especially in the tail region. 

For T=574 K we obtain that the overestimation of 
^MC (*) with respect to Kuu(t) increases, so this be- 
havior suggests that the expression used for Kuc{t) in 
equation ( |22] ) becomes less adequate. Although the ex- 
pression seems acceptable near the triple point, it be- 
comes oversimplified at increasing temperatures and it 
seems necessary to take into account the coupling to 
other modes (i.e., the modes associated with-, the cur- 
rents). In fact, as shown by Shimojo el a£il in their 
MD study for liquid Na, the combined effect of the cou- 
pling terms associated with the currents, gives rise to a 
small and negative contribution to the i^Mc(0 which, if 
taken into account, would lead to a better agreement be- 
tween the K^P\t) and KMu(t)- On the other hand, 
it is also observed that the theoretical i^Mc(^) obtained 
through the present self-consistent scheme follows nicely 
the iifwQ (t), whereas the deficiencies in the viscoelastic 
K^(t) become more marked as the temperature is in- 
creased, with its minimum at t « 0.2 ps becoming deeper 
and taking on negative values. 

The results obtained for the self-diffusion coefficients, 
D, are shown in Table [ij for the three thermodynamic 
states considered in this work. Specifically, the MD re- 
sults have been deduced from the slope of the corre- 
sponding mean square displacement whereas the theoret- 
ical ones where obtained from equation (p3|). First, we 
note that the MD results compare rather well with thp 
different experimental data (INS dataoQ, tracer dataEj) 
for the three temperatures. The theoretical results show 
good agreement with the MD ones for T=470 K whereas 



those for T=526 K and T=574 K are slightly underesti- 
mated. In fact, according to the relation between D and 
the VACF total memory function, see equation (|2^), this 
underestimation of D is a consequence of the above men- 
tioned overestimation of the mode-coupling component 
of K(t) when the temperature is increased. 

In Table [n] we have also included the values for D 
obtained within the viscoelastic approximation. Accord- 
ing to the shortcomings already mentioned for the cor- 
responding memory functions, K v i sc (t), obtained within 
this approximation, when equation ( p3| ) is applied to 
compute the corresponding D, it leads to values which 
overestimate both the MD and experimental results, ex- 
cepting that for T=574 K. Nevertheless, the reasonable 
agreement achieved at the higher temperature is, basi- 
cally due to the cancellation between the overestimation 
and the underestimation of the K v i sc (t) at f k 0.13 ps. 
and t w 0.20 ps. respectively. 

Finally, we end up by noting that the gaussian ansatz 
chosen for the binary part of the K(t), provides an ac- 
curate description of its short time behaviour. We have 
performed calculations by using the squared hyperbolic 
secant in equation ([Hj), and the obtained, selfconsistent, 
results overestimate both the binary and mode-coupling 
parts of the total memory function; in fact it predicts, 
for T=470 K, a self-diffusion coefficient, D=0.55 A 2 /ps 
which is significantly smaller than the corresponding MD 
result. 



C. Dynamic structure factors. 

In this work, both K{t) and F(k, t) are obtained 
within a self-consistent scheme, whereas the self inter- 
mediate scattering function, F s (k, t), is computed within 
the gaussian approximation, see equation (|l^), by us- 
ing the VACF deduced from the self-consistent process. 
Now, by Fourier tranforming F(k,t) and F s (k,t) we get 
S(k,ui) and S s (k,oj) respectively, and by equation ([I]), 
the Stot(k,Lu) is finally obtained. 

In figure |6| we show the dynamic structure factors 
S(k, uj) obtained-from the present formalism, along with 
the MD j-csultsEJ, for T=470 K, and the experimental 
IXS datau, measured at T=488 K, for three wave- vectors. 
The viscoelastic S v i sc (k, uj) are also included for compar- 
ison. The good agreement .between experiment and MD 
has already been remarked^, and supports the adequacy 
of the potential used in order to describe the effective 
interaction among the ions in the liquid. It can be ob- 
served that the overall shape of S(k,uj) is qualitatively 
reproduced by the present theoretical approach, in con- 
trast with the viscoelastic model, which gives correctly 
the peak positions, but fails to describe the overall in- 
dependence of the dynamic structure factor. The behav- 
ior of S(k, uj) is of course a consequence of the time depen- 
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dence of the intermediate scattering functions, F(k,t). 
For instance, the fact that the amplitude of the decaying 
tail of F(k,t) is well reproduced implies that for small 
lo the structure factor will increase significantly, as ob- 
served in the figure. On the other hand, the deficiencies 
in the description of the oscillations of F(k, t) around the 
decaying tail are reflected in the less accurate represen- 
tation of the Brillouin peak of S(k, lo). 

The theoretical results for S to t{k,Lo) as a function of 
lo, at some fixed A-values, for liquid lithium at the three 
temperatures, are shown in figures |^-^|, where they are 
compared with, both the corresponding MD and experi- 
mental dataErH. First, it is observed that the MD results 
show a good agreement with the experimental data, ex- 
cept for some discrepancies appearing at the frequency 
region close to lo « 0. Although these discrepancies could 
be ascribed to the NPA interatomic pair potential used 
in the simulations, this may not be the ultimate reason; 
in fact similac, problems have also been encountered by 
Torcini et a£3 in a similar MD study carried out with 
the interatomic pair potentials proposed by Price et aO. 
Moreover, it must be reminded that an accurate experi- 
mental determination of the spectral features near u) = 
is not a trivial task and, therefore, the above mentioned 
discrepancies between MD and experiment must be taken 
with some caution. 

On the other hand, the comparison MD/theory shows 
that the present theory accounts rather well for the sim- 
ulation results, except in the frecuency region close to 
lo = 0. However, this problem has to be discussed in 
terms of both S s (k,Lo) and S(k,Lo) whose combination, 
according to equation (|l|), gives rise to the Stot(fc,^)- It 
is observed that for k < k p , the main contribution to the 
Stot(k,Lo) comes from the S s {k,Lo), for k w k p it is the 
S(k,Lo) the one which dominates the S , to t(fc,w) whereas 
for k > k p both contributions to the Stot(k, lo) are rather 
similar. Now, when the theoretical S s (k, lo) are compared 
with their MD counterparts, it is found that the theory 
can reproduce the MD results, except at the region close 
to lo = where the theoretical values are always sys- 
tematically lower; in fact, this is a well known deficiency 
of the gaussian approximation used to calculate the self- 
intermediate scattering function. On the other hand, the 
theoretical S(k,co) somewhat overestimates the MD val- 
ues near lo — 0, and therefore, this explains why the 
theoretical St t(k, lo — > 0) are sistematically smaller than 
their MD counterparts except for those fc-values around 
kp. 

In the figures we have also included the theoretical 
results obtained within the simple viscoelastic approxi- 
mation. It is observed that the agreement with the MD 
results is very similar to that obtained with the present 
theory, although now the values for lo — > tend to be 
slightly smaller. This is again, basically, due to the corre- 
sponding S s (k, lo) which when calculated by the gaussian 



approximation incorporates a VACF which gives a bigger 
diffusion coefficient. This implies that F s (k,t) decreases 
faster and leads to a smaller S s (k,LO — > 0). Moreover, 
the values of iSvisc^,^ - ► 0) are also too low, which re- 
inforces the tendency to produce smaller St o t(k,L0 — > 0) 



D. Shear viscosity. 

Figure | shows, for T=470 and 574 K, the MD re- 
sults obtained for the SACF along with its three contri- 
butions. Note that, according to these results, the con- 
tributions from both the potential- kinetic, rjkpft), and 
kinetic-kinetic, rjkk(t), parts of the SACF are negligible 
(less that 5% of the potential-potential part) and there- 
fore it justifies to restrict the theoretical calculations to 
the purely potential part of the SACF, as su ggest ed in 



the theoretical framework presented in section II D 



In the figure, we have also plotted the theoretical r](t) 
along with its binary and mode-coupling components. 
Whereas the binary part alone is quite unable to describe 
the behaviour of ?y(t), the inclusion of the mode-coupling 
component leads to a good overall agreement with the 
simulation data. However, it must be stressed that the 
present theoretical results slightly overestimate the short 
time (t < 0.1 ps) behaviour of i](t) and this is mainly 
due to the binary component; more explicitly, it comes 
from the overestimation of the values for the T n . Again, 
we note that the choice of the squared hyperbolic secant 
ansatz would lead to a worse description of the binary 
part. 

On the other hand, the behaviour for t > 0.1, which 
is completely determined by the mode-coupling compo- 
nent, is rather well described. The values obtained for the 
shear viscosity coefficient, 77, are presented in Table 111. 
First, we stress that the MD results comparej-cather well 
with the corresponding experimental valuescil whereas, 
as expected from figure |^, the theoretical values overesti- 
mate the MD ones because of their short time behaviour. 

Finally, we have also plotted in figure |9|, the results ob- 
tained by the viscoelastic approximation. It is observed 
that it overestimates the short time behaviour whereas 
the long-time behaviour is slightly underestimated; both 
effects counteract so as to produce a shear viscosity co- 
efficient in good agreement with the MD results. 



V. CONCLUSIONS. 

In this paper we have evaluated several dynamical 
properties of liquid lithium at three different thermody- 
namic states. The calculations have been performed both 
by a theoretical framework and by MD simulations, using 
an interatomic pair potential derived within the Neutral 
Pseudoatom method. 
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The MD results have been compared with the exper- 
imental dynamic structure factor, S(q,v), obtained by 
IXSeI, the total dynamic structure factor, Stot(k,uj), ob- 
tained by INSLna, as well as with the experimental data 
for the self-diffusion and shear viscosity coefficients. The 
overall quality of the agreement is good, supporting the 
idea that the interatomic pair potencial derived from the 
NPA method can be used to describe the dynamical prop- 
erties of liquid lithium at thermodynamic states near the 
triple point. 

However, although the MD simulations can always be 
used to evaluate the dynamical properties of a system, 
once the interatomic pair potencial is given, it is also 
important to have a purely theoretical formalism which, 
besides of complementing the MD calculations, will allow 
the analysis of the underlying mechanisms of motion of 
the particles of the system. This is a relevant aspect of 
the present paper, where we have presented a whole theo- 
retical framework which, by incorporating mode-coupling 
concepts, allows the self-consistent determination of sin- 
gle particle properties, as represented by the velocity 
autocorrelation function, its memory function, the self- 
diffusion coefficient and the self-dynamic structure factor, 
as well as collective properties such as the intermediate 
scattering function, the dynamic structure factor, the au- 
tocorrelation function of the non-diagonal elements of the 
stress tensor and the shear viscosity coefficient. The ap- 
plication of this theoretical framework to the study of 
liquid lithium at thermodynamic conditions not far from 
the triple point, leads to results in overall good agree- 
ment with both the corresponding MD results and the 
experimental data. This agreement is particularly sat- 
isfactory at the melting point whereas some deviations 
appear at the highest temperature investigated. 

We have shown that the deficiencies in the VACF mem- 
ory function, observed in previous works and related to 
the use of the viscoelastic model for the intermediate 
scattering functions, are now significantly reduced, lead- 
ing to a much better description of the long time tail of 
K(t). This is basically due to the improved description of 
the intermediate scattering functions, F(k.t), in partic- 
ular, of the overall amplitude of its decaying tail. How- 
ever, although the present theoretical scheme accounts 
for these features of both the single particle and the col- 
lective dynamics, there is still room for improvement in 
some dynamical magnitudes, like F(k,t), for which we 
find some discrepancies in the amplitude of the oscilla- 
tions in the tail as compared with the corresponding MD 
results, which are also reflected in the less accurate de- 
scription of the side peak of the dynamic structure factor. 

Within the present theoretical framework, a key role is 
played by the memory functions; therefore, it is impor- 
tant to provide an accurate description of these magni- 
tudes. Concerning the memory functions given in equa- 
tions (0) and (Fill) , this work shows that its short time be- 



haviour is reasonably described by means of the gaussian 
ansatz. We have verified that, within the approximations 
made in the present work, the hyperbolic secant squared 
ansatz leads to results which are too wide as compared 
with the MD data. 

As a final comment, we note that some authorsE3 have 
used a simplified expresion for the relaxation time of the 
binary part of the second order memory function of the 
intermediate scattering functions. Instead of using the 
Ti(k) which gives the correct sixth moment of S(k,ui), 
they proposed to use the simpler r s (fc), which gives the 
correct sixth moment of S s (k,u>). Certainly, the values 
of Ti(k) oscillate around those of r s (fc)E3, but we have 
found that the influence of using r s (fc) is not negligible, 
giving a noticeably worse description of the intermediate 
scattering functions. 

We end up by signaling the limitations of the formal- 
ism presented here. First, its density/temperature range 
of applicability lies within the region where the relevant 
slow relaxation channel is provided by the coupling to 
density fluctuations, and this ceases to be valid for den- 
sities smaller than those typical of the melting point. For 
these densities, coupling to other modes, like longitudi- 
nal and/or transverse currents, becomes increasingly im- 
portant and we believe that this is the main reason for 
the small deviations observed in the VACF total memory 
function at the two higher temperatures studied. More- 
over, we consider that an improvement in the description 
of the intermediate scattering function, F(k,t), at all 
temperatures, would also require the inclusion of other 
modes in the corresponding second-order memory func- 
tion. 

Second, the theory treats in a different way the self 
intermediate scattering function on one hand and the ve- 
locity autocorrelation function and the intermediate scat- 
tering function on the other, since the mode-coupling ef- 
fects are included for the latters but not for the first. It 
would clearly be desirable a more symmetric treatment, 
and further work is currently performed in that direction, 
and the results will be reported in due time. 
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FIG. 1. Normalized second-order memory function, T(k, t), 
of the intermediate scattering functions, F(k,t), at several 
fc-values, for liquid lithium at T = 470 K. Open circles: 
MD results. Continuous line: present theory. Dash-dotted 
line: binary part, TB(k,t). Dotted line: mode-coupling part, 
TMc(k,t). Dashed line: viscoelastic model. The inset shows 
the normalized first-order memory function as obtained by 
MD (open circles), the viscoelastic model (dashed line) and 
the present theory (continuous line) . 
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FIG. 3. Normalized intermediate scattering functions, 
F(k,t), at several fc-values, for liquid lithium at T = 470 K. 
Open circles: MD results. Continuous line: present theory. 
Dashed line: viscoelastic model. 
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FIG. 2. Same as the previous figure but for T = 574 K. 
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Figure 4. "Dynamic properties of liquid lithium above melting". J. Casas et al 



FIG. 4. Same as the previous figure but for T = 574 K. 
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Figure 5. "Dynamic properties of liquid lithium above melting". J. Casas et a 



FIG. 5. Normalized VACF total memory functions for liq- 
uid Li at three temperatures. Open circles: MD results. Con- 
tinuous line: theoretical results. Dotted line: theoretical re- 
sults for the mode-coupling contribution. Dashed line: vis- 
coelastic approximation. The insets provide a closer compar- 
ison among the mode-coupling contributions from the present 
theory, KMc{t), (continuous line), and from the viscoelastic 
model, #mc(*)j (dashed line). The insets for T=470 K and 
574 K also include the MD results for the mode-coupling con- 
tribution Ky.Q (filled triangles) evaluated by using equation 
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FIG. 6. Dynamic structure factors of liquid Li at T=470 K. 
Continuous line: theoretical results. Dashed line: viscoelastic 
approximation. Open circles: MD, results. Full circles with 
error bars: experimental IXS date 

at T=488 K 
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FIG. 7. Total dynamic structure factors at several fc-values 
corresponding to liquid Li at T= 470 K. Continuous line: 
theoretical results. Dashed line: viscoelastic approximation. 
Open circles: MfjUpesults. Full circles with error bars: exper- 
imental INS datatM 



FIG. 9. Normalised potential part of the stress autocorre- 
lation function, ri(t), for liquid lithium at T—470 and 574 
K. Open circles: MD results. Continuous line: theoreti- 
cal results. Dash-dotted line: binary part. Dotted line: 
mode-coupling component. Dashed line: viscoelastic model. 
The inset shows the MD results for r/(t) (open circles), rj pp (t) 
(continuous line), 10 Xrikp(t) (dotted line) and 10 xr/kk(t) 
(dashed line). 
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FIG. 8. As the previous figure, but for Li at T=526 K and 
574 K. 
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TABLE I. Thermodynamic states studied in this work. 



T (K) 


470 


526 


574 


P (A-') 


0.0445 


0.0441 


0.0438 



TABLE II. Self-diffusion coefficient (in A 2 /ps units), of 
liquid lithium at the thermodynamic states studied in this 
work. D t h, D v [ sc and Dmd are the theoretical, viscoelastic 
and Molecular Dynamics results obtained in this work. 



T(K) 


470 


526 


574 


Ah 


0.689 


0.875 


1.03 


-Dvisc 


0.826 


0.992 


1.14 




0.69 


0.91 


1.11 




0.64±0.04 a 


0.90±0.06 a 


1.08±0.03 



0.69±0.09 b 0.95±0.12 b 1.19±0.14 b 



0.67±0.06 c 0.93±0.07 c 1.16±0.09 c 

a Ref. R 
b Ref. || 
c Ref. |lf 

TABLE III. Shear viscosity (in GPa ps) of liquid lithium at 
the thermodynamic states studied in this work. r; t h and ijmd 
are the theoretical and Molecular Dynamics results obtained 
in this work. 



T(K) 


470 


526 


574 


Vth 


0.589 


0.514 


0.461 


TJvisc 


0.525 


0.481 


0.440 


?)MD 


0.55 


0.48 


0.42 




0.57±0.03 a 


0.50±0.03 a 


0.45±0.03 a 



a Ref. |l] 
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